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Abstract 

We introduce a generalization of the well-known random sequential addition (RSA) process for 
hard spheres in d-dimensional Euclidean space W^. We show that all of the n-particle correlation 
functions of this nonequilibrium model, in a certain limit called the "ghost" RSA packing, can 
be obtained analytically for all allowable densities and in any dimension. This represents the 
first exactly solvable disordered sphere-packing model in arbitrary dimension. The fact that the 
maximal density (l){oo) = 1/2"^ of the ghost RSA packing implies that there may be disordered 
sphere packings in sufficiently high d whose density exceeds Minkowski's lower bound for Bravais 
lattices, the dominant asymptotic term of which is 1/2'^. Indeed, we report on a conjectural lower 
bound on the density whose asymptotic behavior is controlled by 2~('^ ''^^^^'")'^, thus providing 
the putative exponential improvement on Minkowski's 100-year-old bound. Our results suggest 
that the densest packings in sufficiently high dimensions may be disordered rather than periodic, 
implying the existence of disordered classical ground states for some continuous potentials. 
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I. INTRODUCTION 



A collection of congruent spheres in d-dimensional Euclidean space M*^ is called a sphere 
packing if no two spheres overlap. The packing density or simply density of a sphere 
packing is the fraction of space Mf^ covered by the spheres. Hard-sphere packings have been 
used to model a variety of systems, including liquids amorphous and granular media 
P], and crystals [3]. Nonetheless, there is great interest in understanding sphere packings 
in high dimensions in various fields. For example, it is known that the optimal way of 
sending digital signals over noisy channels correspond to the densest sphere packing in a 
high dimensional space jjj. These "error-correcting" codes underlie a variety of systems in 
digital communications and storage, including compact disks, cell phones and the Internet. 
Physicists have studied hard-sphere packings in high dimensions to gain insight into ground 
and glassy states of matter as well as phase behavior in lower dimensions jsf 0]. The 
determination of the densest packings in arbitrary dimension is a problem of long-standing 
interest in discrete geometry ^. 

It is instructive to note that upper and lower bounds on the maximal density 

0max = sup 0(P) (l) 
PCS'* 

exist in all dimensions |^, where the supremum is taken over all packings P in R'^. For 



example, Minkowski [3| proved that the maximal density 0max among all Bravais lattice 



L 



packings for d > 2 satisfies the lower bound 

0max — 2d-l 



^ >m, (2) 



where ({d) = YlT=i ^"'^ Riemann zeta function. One observes that for large values 

of d, the asymptotic behavior of the nonconstructive Minkowski lower bound is controlled 
by 2~'^. Interestingly, the density of a saturated packing of congruent spheres in R'' for all d 
satisfies 

<^ > Yd- (3) 
A saturated packing of congruent spheres of unit diameter and density in JJ'^ has the 
property that each point in space lies within a unit distance from the center of some sphere. 
Thus, a covering of the space is achieved if each center is encompassed by a sphere of unit 
radius and the density of this covering is 2^^^ > 1, which proves the so-called greedy lower 
bound Note that it has the same dominant exponential term as 



A statistically homogeneous (i.e., translationally invariant) packing is completely con- 
figurationally characterized by specifying all of the n-particle correlation functions. For 
such packings in M*^, these correlation functions are defined so that p"'(?„(ri, r2, . . . , r„) is 
proportional to the probability density for simultaneously finding n particles at locations 
ri, r2, . . . , r„ within the system, where p is the number density. Thus, each Qn approaches 
unity when all particle positions become widely separated within M'^, indicating no spatial 
correlations. To date, an exact determination of all of the ra-particle correlation functions 
for a packing has only been possible for d = I'm the special case of an equilibrium ensemble 
of such particles Qj. Observe that in the limit c? — oo, it is known that the pressure of 
an equilibrium hard-sphere fluid is exactly given by the low-density expansion up to the 
second- virial level for a positive range of densities j5|], which implies a simplified form for all 
of the correlation functions j9|. 

We present in Section II a generalization of the well-known random sequential addition 
(RSA) process of hard particles . In a particular limit of this nonequilibrium model 

that we call the "ghost" RSA process, we are able to obtain the Qn for all allowable densities 
exactly for any n and dimension d. The key geometric quantity that determines Qn is the 
union volume of n overlapping exclusion spheres of radius equal to the sphere diameter. 
We show that this construction of a disordered but unsaturated packing realizes the greedy 
lower bound Q. This implies that there may be disordered sphere packings in sufficiently 
high d whose density exceeds Minkowski's lower bound (j2|). Indeed, in Section III, we re- 
port on a conjectural lower bound on the density whose asymptotic behavior is controlled 
by 2~(°-'^'^^^^- -)'^, thus providing the putative exponential improvement on Minkowski's 100- 
year-old bound. Our results lead to the counterintuitive possibility that optimal packings 
in sufficiently high dimensions may be disordered and thus have implications for our funda- 
mental understanding of classical ground states of matter. 



II. GENERALIZED RANDOM SEQUENTIAL ADDITION MODEL 

We introduce a disordered sphere-packing model in that is a subset of the Poisson 
point process and is a generalization of the standard random RSA process. The centers of 
"test" spheres of unit diameter arrive continually throughout M°' during time t > according 
to a translationally invariant Poisson process of density rj per unit time, i.e., 77 is the number 
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of points per unit volume and time. Therefore, the expected number of centers in a region 
of volume fl during time t is rjQt and the probability that this region is empty of centers 
is exp(— ?7f2t). However, this Poisson distribution of test spheres is not a packing because 
the spheres can overlap. To create a packing from this point process, one must remove test 
spheres such that no sphere center can lie within a spherical region of unit radius from any 
sphere center. Without loss of generality, we will set rj = 1. 

There is a variety of ways of achieving this "thinning" process such that the subset 
of points correspond to a sphere packing. One obvious rule is to retain a test sphere at 
time t only if it does not overlap a sphere that was successfully added to the packing 
at an earlier time. This criterion defines the standard RSA process in R"^ 0, which 
generates a homogeneous and isotropic sphere packing in R'^ with a time- dependent density 
4>{t). In the limit t ^ oo, the RSA process corresponds to a saturated packing with a 
maximal or saturation density 0s(oo) = lim(^oo0(t). In one dimension, the RSA process is 
commonly known as the "car parking problem" , which Rehyi showed has a saturation density 
<f)s{oo) = 0.7476 . . . Q|. For 2 < d < oo, an exact determination of 0s(oo) is not possible, 
but estimates for it have been obtained via computer experiments for low dimensions 0]. 

Another thinning criterion retains a test sphere centered at position r at time t if no 
other test sphere is within a unit radial distance from r for the time interval nt prior to t, 
where k is a positive constant in the interval [0,1]. This packing is a subset of the RSA 
packing, and hence we refer to it as the generalized RSA process. Note that when k = 0, 
the standard RSA process is recovered, and when k = 1, a model due to Matern is 
recovered jl^]- The latter is amenable to exact analysis and is the main focus of this paper. 
For any < k < 1, the generalized RSA process is always an unsaturated packing. Figure 
illustrates the differences between the generalized RSA process at the two extremes of k = 
and K = 1. In remainder of this section, we will focus on the 

The time-dependent density 0(t) in the case of the generalized RSA process with k, = 1 
is easily obtained. In this packing, a test sphere at time t is retained only if does not overlap 
an existing sphere in the packing as well as any previously rejected test sphere, which we 
will call "ghost" spheres. The model itself will be referred to as the ghost RSA process. 
An overlap cannot occur if a test sphere is outside a unit radius of any successfully added 
sphere or ghost sphere. Because of the underlying Poisson process, the probability that a 
trial sphere is retained at time t is given by exp(— t>i (l)t), where vi{l) = 7t'^/^/T{1 + d/2) 
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FIG. 1: The addition of four successfully added particles (in the numerical order indicated) in the 
generalized RSA process at the two extremes of k = (left panel) and k = 1 (right panel). In 
both cases, the rejected particles have dashed boundaries. For the case k = 1, a test sphere cannot 
overlap a ghost sphere. Here 3' represents the second attempt to add a third sphere. 

is the volume of a sphere of unit radius. Therefore, the expected time-dependent number 
density p(t) and packing density (pit) = p(t)fi(l/2) at any time t are given by 

pW= /'exp(-„.(im^ ^-^-P<7<^'" . ^(,), l-exP(;^-.(lW . (4) 

Jo ^i(l) ^ 

In the limit t ^ oo, we therefore have that 

p(oo) = lim p{t) = , 0(oo) = lim 0(t) = ^. (5) 

t^oo "^1(1-) t^oo 2 

Observe that the greedy lower bound on the density is achieved in the infinite-time 
limit for this sequential but unsaturated packing, which was pointed out only recently . 
Although the limiting packing density 0(00) = 1/2"^ is far from optimal in low dimensions, 
it is relatively large in high dimensions, as discussed in our concluding remarks. Obviously, 
for any < k < 1, the maximum (infinite-time) density of the generalized RSA packing is 
bounded from below by 1/2'^ (i.e., the maximum density for n = 1). Henceforth, we write 

Vi = Vi{l). 

The derivation of the expression of g2{r',t) is actually a simple extension of the afore- 
mentioned one for p(t). Two test spheres that arrive at times ti and ^2 and whose centers 
are separated by a distance r can only be retained if no other test spheres arrived before 
ti and t2, respectively (see Fig. |2)). Thus, the key geometrical object is the union volume 
f2(r) of two spheres of unit radius whose centers are separated by a distance r, which can 
be expressed in terms of the intersection volume v^^'^{r) [13] between two such spheres via 
the relation 

V2{r) = 2^1 -<*(r). 



FIG. 2: Left panel: Relevant subvolumes for two overlapping spheres of unit radius associated 
with the arrivals of two test spheres. The labels refer to distinct, nonoverlapping regions. Right 
panel: Relevant subvolumes for three overlapping spheres of unit radius associated with the arrivals 
of three test spheres. 



For r > 2, there is no volume common to two such spheres (f. 



intl 



0) and therefore 



92{j';t) = 1, i.e., pair correlations vanish. However, if 1 < r < 2, the two spheres have a 
common volume and 



p\t)g2{r-t) 



exp - U[vi - vf{r)] - t2[vi - vT{r)] - max (ti, t2)vT{r) 

'0 JQ 

ri /•t2 
2 / dt2 exp[-t2Vi] / dtiexp - ti[vi - v''^\r)] 
Jo Jo 



dtidti 



V2[r) - Vl 



1 



3-to2(r) 



Vl V2[r) J 

In relation the terms within the first three brackets are the distinct volumes of the 
regions labeled 1, 2 and 12 in the left panel of Fig. El Therefore, the time-dependent pair 
correlation function for all r and t is given by 



p\t)g2{r;t) 



2e(r - 1) 

V2{r) - Vl 



1 - e 



Vl 



I _ Q-tv2{r) 

V2{r) 



(7) 



where Q{x) is the unit step function, equal to zero for x < and unity for x > 1. It is useful 
to note that at small times or, equivalently, low densities, formula (j^ yields the asymptotic 
expansion (j){t) = t — + O(t^), which when inverted yields t = + 2'^~^(j? + 

Substitution of this last result into (|7j) gives 



^2(r;0) = e(r-l) + O(03), 



(8) 



which implies that g2{v'-, 0) tends to the unit step function 0(r — 1) as ^ for any d. 

In the limit t —>■ oo, we have from that p'^{oo)g2{r; oo) = 20(r — l)/[fif2(r)] or, using 
(3), 

, 2e(r-l) 



where ^2(1^) = V2{r)/vi. The radial distribution function g2{r; 00) is plotted in Fig. Olfor the 
first five space dimensions. Because /32(r) is equal to 2 for r > 2, g2{r; 00) = 1 for r > 2, i.e., 
spatial correlations vanish identically for all pair distances except those in the small interval 
[0, 2). Even the positive correlations exhibited for 1 < r < 2 are rather weak and decrease 

n 

exponentially fast with increasing dimension |9|, i.e., g2{r; 00) tends to the unit step function 
as d ^ 00, i.e., beyond the hard core (a constrained correlation), spatial correlations vanish. 

Matern originally gave an expression for the time-dependent density (j){t) and and a 
formal expression (as opposed to explicit expression for any d) for the time-dependent radial 
distribution function g2{r;t) when k = 1 using a completely different approach. However, 
he did not consider obtaining any of the higher-order correlation functions. 
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FIG. 3: Radial distribution function for the first five space dimensions at the maximum density 
(p = 1/2'^ for the generalized RSA model with k = 1, i.e., the "ghost RSA process." 

Let us now derive the time-dependent triplet correlation function 5'3(ri2, ri3; 00). Here the 
relevant geometrical object is the union volume f3(ri2, ris, r23) of three spheres of unit radius 
whose centers are separated by the distances ri2, ri3 and r23, which can be expressed in terms 
of the intersection volume v^^'^{ri2, ri3, r23) between three such spheres via the relation 

^3(ri2, ri3, r23) = 3v, - vT\n2) - v^^'in^) - <*(r23) + vl^\n2, ri3, r23). (10) 

Whenever there is no overlap between the three spheres, g^ = 1, i.e., triplet correlations 
vanish. On the other hand, whenever the spheres overlap such that each pair distance 
is greater than or equal to unity, there are triplet correlations. In such situations, it is 
convenient to introduce the time-dependent triplet function 

F{n2, ri3, r23; h, t2, h) = -h[v, - i;r*(ri2) - ^^^^(^3) + vl^\r,2, ri3, r23)] 
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-^2^1 - <*(ri2) - <*(^23) + <*(ri2, ri3, rsa)] - h[vi - v^^ns) - v'^'^'ir^s) + vl^'in^, r^, ras)] 
-max(ti,t2)K*(n2) -<*(ri2,ri3,r23)] -max(ti,t3)K*(n3) -<*(n2,ri3,r23)] 
- max (t2, t3)br*(^23) - vl''\ri2, ri3, r23)] - max (ti, t2, t3)t^r*('"i2, ^"13, '^23)- 



The terms within the first three brackets are the volumes of the regions labeled 1, 2 and 3 
in the right panel of Fig. |21 The terms within the fourth through sixth brackets are the 
volumes labeled 12, 13 and 23 in the right panel of Fig. |21 Of course, the region labeled 123 
denotes the intersection volume of three spheres. The triplet correlation function at time t 
is given by 



P^Wfl'3('^12,n3,'^23;^) 



t ft rt 



exp 



dtidt2dt3 



iQ Jo Jo 

and therefore at infinitely large times we have, using (0 Q and (fT^ . that 



P^{oo)g3{ri2,ri3,r23; 00) 



t3 



dt3exp[-t3Vi] dt2exp - ^2^ - fa" (r23)] 



t2 



X / dh exp - t,[v, - <*(ri2) - vT\n3) + < (ri2, ^3, r23)] 
Jo 

+ 2 dtiexp[-tivi] / rft3exp - t3[vi - vi^'^ns)] 

Jo Jo 

X / rft2exp -t2K -^^2"*('^12) -^^2"*('^23) +<*(ri2,ri3,r23)] 











t2 



+ 2 dt2exp[-t2Vi] / dtiexp -ti[fi-f2 (^2 



ti 



X / dt3 exp 
'0 

2 

flf3('^12,n3,'^23) 



- ^3^1 - <(ri3) - <(r23) + <(ri2, ri3, r23)] 



+ 



1 



+ 



.^2(^2) ' ^2(^^13) ' ^2(^^23)]' 

Combination of (0 © and (fT^ yields the following expression for the triplet correlation 
function for arbitrary positions at infinitely large times: 
0(ri2 - l)e(ri3 - l)e(r23 - 1) 



5'3(ri2,ri3; 00) 



5'2(ri2;oo) + 5(2(ri3;oo) + g2{r23; ooll^) 



f33{ri2,ri3,r23) 

where /53(ri2,ri3,r23) = f 3(^^12, ^3, ^^23)/^! and 5(3(ri2, ri3; 00) = fi'3(ri2, ri3, r23; 00). 

A similar analysis reveals that the four-particle correlation function in the limit t ^ 00 
is given by 



5'4(ri2,ri3,ri4; 00) 



n-<,e(r.,-i) 

/54(ri2,ri3,ri4; 00) 
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fi'3(ri2, ri3; oo)+5(3(ri2, r^; oo)+5(3(ri3, Tm; cx))+5(3(r23, r24; 00) 

(14) 
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By induction, the n-particle correlation function for arbitrary positions at infinitely large 
times is given by 

5f„(ri, . . . , r„; oo) = -— > gn-i[Qu oo) , (15) 

where the sum is over all the n distinguishable ways of choosing n — 1 positions from n posi- 
tions ri, . . . r„ and the arguments of Qn-i are the associated n — 1 positions, which we denote 
by Qi. Moreover, /5„(ri2, ris, . . . , ri„) = t;„(ri2, ri3, . . . , ri„)/t;i, where t;„(ri2, ris, . . . , ri„) 
is the union volume of n congruent spheres of unit radius whose centers are located at 
ri, . . . , Tn, Vij = Tj — Ti for all 1 < z < j < n and Vij = \rij\. 
It can be shown 9] that in the limit d oo and for = 1/2'^ 

n 

5(„(ri2, . . . ,ri„; oo) ~ JJfl'2(?^»i; oo), (16) 

where g2{r; oo) ~ 0(r — 1). We see that unconstrained spatial correlations vanish asymp- 
totically. Specifically, (i) the high-dimensional asymptotic behavior of g2 is the same as the 
asymptotic behavior in the low-density limit for any d [cf. ((HI)], i.e., unconstrained spatial 
correlations, which exist for positive densities at fixed d, vanish asymptotically for pair dis- 
tances beyond the hard-core diameter in the high- dimensional limit; and (ii) Qn for n > 3 
asymptotically can be inferred from a knowledge of only the pair correlation function g2 
and number density p. These two asymptotic properties, which we have called the decor- 
relation principle 9|, apply more generally to any disordered packing, as discussed in Ref. 
j^. Asymptotically, unconstrained correlations vanish (i.e., statistical independence is es- 
tablished) because we know from the Kabatiansky and Levenshtein asymptotic upper bound 
on the maximal density 0max of any sphere packing that the density must go to zero at least 
as fast as 2-^-^^^^'^ for large d Q. 



III. DISCUSSION 

The fact that the maximal density 0(oo) = 1/2'^ of the ghost RSA packing coincides with 
the greedy lower bound (jHI) strongly suggests that there are saturated disordered packings 
;hat have larger densities, i.e., the greedy lower bound is a weak bound for saturated packings 
isl ] . This implies that there may be disordered sphere packings in sufficiently high d whose 
density exceeds Minkowski's lower bound ^ for Bravais lattices, the dominant asymptotic 

9 



term of which is 1/2'^. Our results already give insight into this fascinatingpossibility. 



Fo. example, con.de. ..caUed ^oa^d la.ice in . di^en^o.^. wHich . 

a (i-dimensional generalization of the optimal (densest) face-centered cubic lattice in three 
dimensions, and thought to be the optimal packing for d = A and = 5. Its packing density 
= 7r'^/^/[r(l + (i/2)2^'^+^)/^] exponentially decreases with increasing d (because it quickly 
becomes unsaturated) and falls below the ghost-RSA-process value of 1/2'^ for the first time 
at (i = 28 The ratio of densities of the ghost RSA process to the checkerboard at 

d = 100 is given by 4>ghost/<Pchecker ~ 7.5 x 10^^. Although both packings are unsaturated in 
such high dimensions, the fact that g2{r) for the ghost RSA process is effectively uniform 
(unity) for all r > 1 but for the checkerboard lattice involves Dirac delta functions of weak 
strength at widely spaced discrete distances explains why the former is enormously denser 
than the latter. 

Over the last century, many extensions and generalizations of Minkowski's lower bound 
(j2)) have been obtained ^ , but none of these investigations have been able to improve upon 
the dominant exponential term 2~'^. In another work j^, we will present comprehensive 
rigorous evidence that this exponential improvement may be provided by considering specific 
disordered sphere packings. Here we simply sketch the procedure leading to this putative 
improvement over Minkowski's lower bound. The basic ideas underlying our new approach 
to the derivation of lower bounds on 0max were actually described in our earlier work ^| in 
which we studied so-called (72-invariant processes. A g2-invariant process is one in which a 
given nonnegative pair correlation 5'2(r) function remains invariant for all r over the range 
of densities 

< < 0*. (17) 

The terminal density 0* is the maximum achievable density for the 5f2-invariant process 
subject to satisfaction of certain necessary conditions on the pair correlation. In particular, 
we considered those "test" g2{ry?> that are distributions on R"^ depending only on the radial 
distance r. For any test 02(1^), we want to maximize the corresponding density satisfying 
the following three conditions: 
(i) 92(1^) > for all r, 
(i) g2{r) = for r < 1, 
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S{k) = l + p (27r)2 / /'^/i(r) '^'^^\a]l)J dr > for all k, 

Jo ( fvT j 



(kr) 

where h{r) = (72 (r) — 1 is the total correlation function. Condition (i) is a trivial consequence 
of the fact that (72 is a probability density function. Condition (ii) is just the hard-core 
constraint for spheres of unit diameter. Condition (iii) states that the structure factor S{k) 
in d dimensions must be nonnegative for all k. When there exist sphere packings with g2 
satisfying conditions (i)-(iii) for </> in the interval [0,0*], then we have the lower bound on 
the maximal density given by 

(/"max > 4>*- (18) 

It is rather remarkable that the optimization problem defined above is identical to one 



formulated by Cohn 



. Specifically, it is the dual of the primal infinite-dimensional linear 



program that Cohn employed with Elkies to obtain upper bounds on the maximal 
packing density. Thus, even if there does not exist a sphere packing with g2 satisfying 
conditions (i)-(iii), the terminal density 0* can never exceed the Cohn-Elkies upper bound 
and, more generally, our formulation has implications for upper bounds on 0max- 

In addition, to the structure factor condition, there are generally many other conditions 
that a pair correlation function corresponding to a poi nt p rocess must obey [20]. One such 
additional necessary condition, obtained by Yamada j2l|, is concerned with the variance 
cr^(fi) = {[N{Qy — {N{Q))y), in the number N{Q) of particle centers contained within a 
region or "window" C M'^: 

a^{n) = p\n\ h + P [ h{r)dr\ > 9{1 - 9), (19) 

where 6 is the fractional part of the expected number of points p\Q\ contained in the window. 
This is a consequence of the fact that the number of particles in any window must be an 
integer. 

In Ref. 0. a test fan,,,, of had been cons,de,ed. wMcK .neorp.ated 

the known features of core exclusion, contact pairs, and damped oscillatory short-range order 
beyond contact that are features intended to describe disordered jammed sphere packings 
for d = 3. However, because of the functional complexity of this test g2, the terminal density 
could only be determined numerically. The general optimization procedure outlined above 
was employed in Ref. Q] to obtain analytical estimates of the terminal density in high 
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dimensions that together with the following conjecture provide the putative exponential 
improvement on Minkowski's lower bound on 0max: 

Conjecture 1: A hard-core nonnegative tempered distribution (72 (r) is a pair correlation 
function of a translationally invariant disordered sphere packing in at number density 
p for sufficiently large d if and only if 5'(k) > 0. The maximum achievable density is the 
terminal density 0*. 

In other words, g2{f) that meets the conditions (i) - (iii), at or above a critical dimension 
packings exist with such a (72- A disordered packing in R'' is defined in Ref. 9 to be one 
in which the pair correlation function (72 (i") decays to its long-range value of unity faster 
than |r|^'^~^ for some e > 0." Employing the aforementioned optimization procedure with a 
certain test function (72 and Conjecture 1, we obtain in what follows conjectural lower bounds 
that yield the long-sought asymptotic exponential improvement on Minkowski's bound. An 
important feature of any dense packing is that the particles form contacts with one another. 
Experience with disordered jammed packings in low dimensions reveals that the contact or 
kissing number as well as the density can be substantially increased if there is there is a 
low probability of finding noncontacting particles from a typical particle at radial distances 
just larger than the nearest-neighbor distance. It is desired to idealize this small-distance 
negative correlation (relative to the uncorrelated value of unity) in such a way that it is 
amenable to exact asymptotic analysis. Accordingly, a test radial distribution function was 
considered in Ref. 8 in which there is a gap between the location of a unit step function and 
the delta function at finite i.e., 



where s(r) is the surface area of a d-dimensional sphere of radius r and Z is a parameter, 
which is the average contact or kissing number, and unity is the sphere diameter. The 
expression contains two adjustable parameters, a > 1 and Z, which must obviously be 
constrained to be nonnegative. 

Before reporting the main results of this optimization, it is instructive to examine the 
test function pO|l for two special cases: (1) one in which cr = 1 and Z = and (2) the 
other in which a = 1 and Z > (which were first considered in Ref. 17). In the first 
special instance, there are no parameters to be optimized here, and the terminal density 0* 
is given by 0* = ^. It is simple to show that the Yamada condition ((fTIJ|) is satisfied in 




(20) 
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any dimension for < (f) < 2~'^. We already established in the previous section that there 
exist sphere packings that asymptotically have radial distribution functions given by the 
simple unit step function for < 2^*^. Nonetheless, invoking Conjecture 1 and the obtained 
terminal density, implies the asymptotic lower bound on the maximal density is given by 

0max > ^, (21) 

which provides an alternate derivation of the elementary bound (jH)). Using numerical simu- 
lations with a finite but large number of spheres on the torus, we have been able to construct 
particle configurations in which the radial distribution function is given by the test function 
()20|) with a = 1 and Z = in one, two and three dimensions for densities up to the terminal 



density [22, |23[. The existence of such a discrete approximation to this test g2 is sugges- 
tive that the standard nonnegativity conditions may be sufficient to establish existence in 
this case for densities up to 0*. In the second special case (cr = 1 and Z > 0) and under 
the constraint that the minimum of S{k) occurs at k = 0, then we have the exact results 
0* = and = |, where is the optimized average kissing number. The Yamada 
condition f ()19p is violated here only for d = 1 and becomes less restrictive as the dimension 
increases from d = 2. Interestingly, we have also shown via numerical simulations that there 



jv this test function in 
23l |. This is suggestive 



exist sphere packings possessing radial distribution functions given 
two and three dimensions for densities up to the terminal density 
that the Conjecture 1 for this test function may in fact be stronger than is required. In the 
high-dimensional limit, invoking Conjecture 1 and the obtained terminal density, yields the 
conjectural lower bound 

0max > (22) 

This lower bound provides the same type of linear improvement over Minkowski's lower 
bound as does Ball's rigorous lower bound j3] obtained using a completely different ap- 
proach. 

Now let us consider the problem when both a and Z in fl20j) must be optimized. The 
presence of a gap between the unit step function and delta function will indeed lead asymp- 
totically to substantially higher terminal densities. For sufficiently small d {d < 200), the 
optimization procedure is carried out numerically jQ]. The Yamada condition (fTT?|l is vio- 
lated only for d = 1 for the test function ()20p for the terminal density 0* and associated 
optimized parameters a^, and Z^^ = {2a^(j)^Y — 1. One can again verify directly that the 
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Yamada condition becomes less restrictive as the dimension increases from d = 2. How- 
ever, although the test function ()20|) for d = 2 with optimized parameters 0* = 0.74803, 
(T^, = 1.2946 and Z^, = 4.0148 satisfies the Yamada condition, it cannot correspond to a 
sphere packing because it violates local geometric constraints specified by a* and [9]. To 
our knowledge, this is the first example of a test radial distribution function that satisfies 
the two standard non- negativity conditions (i) and (iii) and the Yamada condition ()19|). but 
cannot correspond to a point process. Thus, there is at least one previously unarticulated 
necessary condition that has been violated in the low dimension d = 2. As is the case with 
the Yamada condition (fT!I|) . this additional necessary condition appears to lose relevance in 
low dimensions because we have shown that there is no analogous local geometric constraint 
violation for d > 3. For d < 56, the terminal density lies below the density of the densest 
known packing (a Bravais lattice) However, for d > 56, can be larger than the density 
of the densest known arrangements, which are ordered. Our numerical results for d between 
3 and 200, reveal exponential improvement of the terminal density 0* over the one for the 
gapless case, where 0* = (rf + 2)/2'^+-'^. 

For large d, an exact (but nontrivial) asymptotic analysis can be performed 0], yield- 
ing the optimal terminal density. This result in conjunction with Conjecture 1 yields the 
conjectural asymptotic lower bound 

3.276100896d^/6 3.276100896 rf^/^ 



2[3-log2(e)]d/2 20-7786524795... d 



(23) 



This putatively provides the long-sought exponential improvement on Minkowski's lower 
bound. We call this a conjectural lower bound because it relies on Conjecture 1 being true, 
which a number of results support. First, the decorrelation principle states that uncon- 
strained correlations in disordered sphere packings vanish asymptotically in high dimensions 
and that the Qn for any n > 3 can be inferred entirely from a knowledge of p and g2- Second, 
the necessary Yamada condition appears to only have relevance in very low dimensions. 
Third, we have demonstrated that other new necessary conditions also seem to be germane 
only in very low dimensions. Fourth, we recover the form of known rigorous bounds [cf. 
(21) and (22)] in special cases of the test radial distribution function (20) when we invoke 
Conjecture 1. Finally, in these two instances, configurations of disordered sphere packings 
on the torus have been numerically constructed with such g2 in low dimensions for densities 
up to the terminal density |2^. I23I. 
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A byproduct of the bound ()23|) is the conjectural asymptotic lower bound on the maximal 
kissing number Q 

^max > ~ 40.24850787 (f^^ 2[i°S2(e)-i]d/2 _ 40.24850787 d^^^ 2°-22i3475205... ^^4) 

This result is superior to the best known asymptotic lower bound on the maximal kissing 
number of 20-2075...rf 

The work described above suggests that the densest packings in sufficiently high dimen- 
sions may be disordered rather than periodic, implying the existence of disordered classical 
ground states for some continuous potentials. In fact, there is no fundamental reason why 
disordered ground states are prohibited in low dimensions 26]. A case in point are the 
"pinwheel" tilings of the plane, which possess both statistical translational and rotational 
invariance [27 1. 

This work was supported by the National Science Foundation under Grant No. DMS- 
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